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ABSTRACT 

A homogeneous sample of ~2200 low redshift disk galaxies with both high sensitivity long-slit optical 
spectroscopy and detailed /-band photometry is used to construct average, or template, rotation curves 
in separate luminosity classes, spanning 6 magnitudes in /-band luminosity. The template rotation 
curves are expressed as functions both of exponential disk scale lengths and of optical radii Ropt, 
and extend out to 4.5-6.5 tr, depending on the luminosity bin. The two parameterizations yield 
slightly different results beyond Ropt because galaxies whose Ha emission can be traced to larger 
extents in the disks are typically of higher optical surface brightness and are characterized by larger 
values of //opt/^’d- By either parameterization, these template rotation curves show no convincing 
evidence of velocity decline within the spatial scales over which they are sampled, even in the case 
of the most luminous systems. In contrast to some previous expectations, the fastest rotators (most 
luminous galaxies) have, on average, rotation curves that are flat or mildly rising beyond the optical 
radius, implying that the dark matter halo makes an important contribution to the kinematics also 
in these systems. The template rotation curves and the derived functional fits provide quantitative 
constraints for studies of the structure and evolution of disk galaxies, which aim at reproducing the 
internal kinematics properties of disks at the present cosmological epoch. 

Subject headings: galaxies: kinematics and dynamics — galaxies: spiral — galaxies: structure — 
cosmology: observations — dark matter 


1. INTRODUCTION 

The determination of the rotational characteristics 
of galaxy disks is fundamental to understanding the 
role that dynamics plays in galaxy formation and evo¬ 
lution over cosmic time. Over the last decade, a 
large number of studies have addressed the issue of 
disk formation and evolution in current hierarchical 
cosmologies by means of analytical or semi-analytical 
theoretical models fe.g..|D alcanton. Soergel, Sl Su mmqrj 
199^jMo^MaOj&Wlii^l99^jMao^^Oj_^^Vhit^l993 

a lFirniaiTTX: Avila,-Rees3l2flflfriTOTrTieT7^oschll2flflir 

and num erical smooth particle hydrodynamic sim¬ 
ulations Te.g-.. Iv an den Bosch et aT ] l2nnl lAbadi et al.1 
120031: iGovernato et al. I 1 200411 . Other groups have 
investigated more specifically the characteristics of 
the density profiles of dark matter halos, which 
directly determine the observable rotation curves 
(RGs ) of present-day disks, using iV-bodv simulations 

('e.g..lNavarro. Frenk. fc WhitelI99dll997l:lBullock et all 

l2001alb[l . A detailed characterization of the kinematic 

properties of observed disk galaxies at low redshift is of 
key importance to validate such models and place con¬ 
straints on the properties of dark matter halos. In this 
context, special attention has been recently devoted to 
dwarf and low surface brightness galaxies, i.e. systems 
with an internal kinematics that is thought to be dom¬ 
inated by dark matter even at small radii, thus allow- 
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ing to b est probe the inner structure of their p arent ha¬ 
los ('e.g.. lSDekkens. Giovanelli. fc Havne J2005I and refer¬ 
ences therein). For high surface brightness galaxies, the 
kinematics is easier to measure, but the presence of a 
more prominent stellar disk complicates the analysis of 
the RG, whose decomposition into a dark matter and a 
disk component is generally not unique (see, for instance, 
l.limenez. Verde, fc Ohll2003t iDiitton et al.ll200M and ref¬ 
erences therein). 

The derivation of average RGs representative of the 
kinematic properties of a large number of local spirals is 
also important on its own, to provide a standard refer¬ 
ence against which similar, distant samples can be com¬ 
pared. Furthermore, statistical descriptors of the am¬ 
plitude and radial variation of rotation velocity within 
disks would benefit all the applications that rely on the 
assumption of a RG model, such as determinations of ro¬ 
tational widths from spectroscopic data that lack velocity 
or spatial resolution. Examples of the latter include in¬ 
ternal kinematic studies of intermediate redshift spirals 
through spatially resolved optical spectroscopy, where 
synthetic RGs are needed in order to estimate the full 
velocity width (e.g. , rtet_^ ]ll997t |Simard fc Pritched 
Il999l[^eg1er et 

Beginning perhaps with lRob'ert^TRota l|1973ll . many 
previous works have explored the variation in rotation 
curve form among galaxies. As larger samples of rota¬ 
tion curves have been amassed, attempts have been made 
to produce more quantitative descriptions of empirically- 
derived average, or template RGs. Notably, average RGs 
binned by luminosity and with radial distances expressed 
in units of the optical radius /?opt defined as the radius 
encompassing 83% o f the total integrated light wer e pre¬ 
viously presented by iPersic. Salucci. fc Stell lll996L here¬ 
after PSS96) for a sample of 616 spiral galaxies. This 
definition of Ropt is equivalent to 3.2 tr for a pure ex- 
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ponential disk. PSS96 showed that their average RCs 
were well described by an analytical form, obtained as 
the combination in quadrature of a dark matter halo 
and an exponential disk, with fit coefficients depend¬ 
ing on galaxy luminosity only. Based on these results 
they claimed the existence of a universal relation, the 
universal rotation curve (URC), whereby the shape and 
amplitude of an observed RC at any radius is com¬ 
pletely determine d by the galaxy lumin osity, as origi¬ 
nally proposed by ^erai^^^alu^ According to 

the URC predictions, the most luminous systems have 
RCs that peak at ~0.8 i?opt and decline beyond that 
radius, implying a minor contribution of the dark mat¬ 
ter halos to the internal kinematics of those galaxies. It 
should be also noticed, however, that PSS96’s average 
RC tracings beyond i?opt are based on linear extrapo¬ 
lations of outer velocity gradients as functions of lumi¬ 
nosity mostly obtained from H I data. Several studies 
have discus sed the inadequacy of the URC parameteri¬ 
zation (e.g., jCourteau|li)97 yyert^eiJl997tlWillickll999l: 
iGarrido. MarceliiTTrAmrainll^OO^ 

The purpose of this paper is to provide a set of tem¬ 
plate RCs in bins covering a wide range of galaxy lumi¬ 
nosity and based on an extensive set of high sensitivity 
long-slit optical RCs, for which a homogeneous body of 
high quality /-band photometry is also available. Our in¬ 
tent is not to propose a new universal relation, but rather 
to obtain a more reliable characterization of the average 
RCs of late type, high brightness disk galaxies by taking 
advantage of a much larger data set, thereby negating the 
need for extrapolation within the optical disk. In order 
to allow a direct comparison with PSS96’s average RCs 
and URC models, we provide template RCs parameter¬ 
ized as functions of //opt- Furthermore, we also present 
a set of templates which parameterize the distance from 
the center of a galaxy by means of the exponential disk 
scale length rd, a quantity naturally adopted in models 
of galaxy disks. 

This work makes use of a large spiral galaxy photomet¬ 
ric and spectroscopic dataset, dubbed SFU—I- (see §EI), 
compiled to investigate the internal kinematics of disk 
galaxies at low redshifts, as well as to study the peculiar 
velocity field in t he local universe via a pplication of the 
Tully-Fisher ('TF MTuIIv fc Fish'^ll977ll method. Details 
on the extraction of related quantities from long-slit spec¬ 
tra, including our RC folding technique and algorithm for 
the measurement of rot ational velocities for TF applica¬ 
tions, can be found in iCatinella. Havnes. fc Giova,ne1lil 
lf2fm hereafter CHG05). 

This paper is organized as follows. In §|21we describe 
our data sample and illustrate the procedure used to de¬ 
rive template RCs. The results obtained using the pa- 
rameterizations in terms of exponential disk scale lengths 
or optical radii are presented and compared in § 13 In 
order to explain systematic differences between the outer 
slopes of the two sets of template RCs, the correlation 
between //opt/?'d ratio and RC extent observed for the 
galaxies in our sample is also investigated in some detail. 
The impact of internal extinction on the inner slopes of 
RCs is analyzed in §0] where template RCs are computed 
for separate intervals of galaxy inclination. The compar¬ 
ison with previously published results, most notably the 
URC models, is discussed in § |3 and our conclusions 
summarized in § ini A value for the Hubble constant of 


Ho = 70 km s ^ Mpc ^ is assumed throughout this work. 

2. DATA SAMPLE AND TEMPLATE RC DERIVATION 

A compilation of /-band photometry and rotational 
width data (obtained from long-slit optical rotation 
curves and/or spatially integrated H I-line profiles) for 
several thousand spiral galaxies has been assembled by 
our group with the principal purpose of mapping the pe¬ 
culiar velocity field in the local Universe. This sample, 
referred to as SFI- |--|-, includes several already published 
datasets ('CHGOSjSnringob et al.l200,^IVogt et al.l2004 

iDale fc Giovanellil 1200(1 and references therein) along 

with recently obtained observations; detai ls of this com¬ 
pilation will be fully discussed elsewhere l|Master^l200,^ 
iSoringobl 1200,^ . Most notably, SFU—I- includes the 
southern spiral galaxy s ample of Mathewson and collab¬ 
orators (MFB sample: [Ma^he wsou. Ford, fc HuchhornI 
Il992t iMathewson fc FordI 1 199(11 ! The /-band isophotal 
photometric fits and RCs derived from long-slit optical 
spectroscopy for the MFB sample are available to us in 
digital form, and were reprocessed to achieve consistency 
with the rest of SFI-I--I- (photometric quantities and rota¬ 
tional widths ha ve been extracted ac cording to the meth¬ 
ods described in iHavnes et al.llT999l and CHG05, respec¬ 
tively) . 

In this work, we use the subset of galaxies in the SFI-|--f 
sample for which high quality long-slit Ha/[N II] RCs 
and /-band photometry are available. Long-slit spectra 
are of high quality when an extended RC can be reli¬ 
ably extracted, i.e. when the Ha emission line has high 
signal-to-noise, is extended, and allows to adequately 
trace the shape of the RC across its spatial extent. Anal¬ 
ogously, /-band images are of high quality when photo¬ 
metric parameters can be reliably extracted (i.e., pho¬ 
tometric quality, no stars superimposed on the galaxy, 
object not close to the edge of the frame, etc. Also dis¬ 
carded are cases in which the images are of high quality 
but the fits are poor because of the presence of a bar, 
asymmetries, or other problems). Template RCs were 
obtained by binning the galaxies in /-band luminosity 
intervals, and computing the average RC in each bin as 
described below. Two separate sets of solutions were de¬ 
rived by adopting different spatial parameterizations, in 
which the radial coordinates of the RCs are expressed in 
units of the /-band disk scale length, rd, or the optical ra¬ 
dius, //opt (i.e., the radius within which 83% of the total 
/-band light is included, as in PSS96. Notice that //opt is 
a cumulative quantity and not an isophotal radius). The 
template RCs were then fitted with the analytic function: 

UpE(r) = Uo(l - e-’'/’-^^)(l + ar/rpE) (1) 

referr ed to as the Polyex model llGiovanelli fc HavnesI 
| 2003 ), where Vq, tpe, and a respectively determine the 
amplitude, the exponential scale of the inner region, and 
the slope of the outer part of the RC (for template RC 
fitting, r and rpE are expressed in units of rd or //opt)- 
Equation ^ is an empirical expression that fits well a 
large variety of RC shapes (including those declining at 
large radii), as discussed in CHG05. Alternative fitting 
functions have been proposed by o t her authors (e.g ., 
lCourteaulll997t iKravtsov et alJir998t iVogt et M]l2004jl : 
since they are all empirical, trade-offs between simplic¬ 
ity and detail result. A comparison between Polyex and 
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other fitting expressions is beyond the scope of this paper 
and will not be discussed here. 

All the individual RCs used for the derivation of 
the template curves were fitted using Equation ^ and, 
although their Polyex parameters were not used in this 
work, a quality code based on the goodness of the fit was 
used to select those included in the final sample. For 
reasons that will be clarified later, the following objects 
were also not used to construct template RCs: (a) 
galaxies with inclinations to the line of sight i > 80° (see 
also ®; (b) galaxies with absolute /-band magnitude 
lying outside the [-24.0,-18.0] interval; (c) RCs with 
spatial extent smaller than 2 rj or 0.6 i?opt) for the rd or 
Ro-pt parameterizations, respectively; and (d) RCs with 
average velocity (computed over the spatial interval Ar, 
as defined in [3] below) incompatible with the value 
obtained using all the RCs in the same luminosity bin 
(over the same Ar). Most of the rejected objects fell 
into category [a] (472 out of 2777 galaxies in the initial 
sample). Also, because of the selection criterion [c], the 
sizes of the final data sets used to derive rd and i?opt 
template RCs, hereafter referred to as ‘Vd and /?opt 
samples”, were slightly different, including 2155 and 
2169 galaxies respectively. 

After binning the galaxies by luminosity, template RCs 
were obtained as follows: 

(1) Folding - Each RC was folded about its center of sym¬ 
metry, as determined by fitting a URC model (PSS96) 
to the unfolded RC. As pointed out elsewhere (see ref¬ 
erences in §^, the URC parameterization is inadequate 
to model individual RCs, but it approximates the over¬ 
all RC shape well enough for the purpose of determining 
its center of symmetry, as argued in CHG05, where our 
folding technique is described in detail. 

(2) Rescaling and Resampling - The radial variable of 

each folded RC was scaled by the rd (i?opt) value of the 
galaxy, and the RC amplitude was deprojected to the 
edge-on view. Inclinations were obtained from elliptici- 
ties measured from isophotal fits to the /-band images, 
and corr ected for seeing effects , following the prescrip¬ 
tions of iCiovanelli et ahl RC amplitudes were 

then resampled on a fixed spatial grid of spacing 0.1 rd 
(0.03 i?opt) using a linear interpolation including the four 
observed points nearest to the grid position. This assures 
that RCs of nearer galaxies, with denser spatial sampling, 
do not bias the average RCs. We have tested that the 
results are not sensitive to the choice of the spatial bin 
centers. Only the resampled RCs were used after this 
point. 

(3) Normalization - Obtaining an average RC for each 
luminosity class requires attention to the fact that the 
RC of each galaxy is sampled out to a different maxi¬ 
mum radius //max- To determine the average amplitude 
of the template RC of a given luminosity class, we thus 
used only the section Ar between 2r/rd and 3r/rd (0.6- 
0.9 r///opt) of each RC. Over this interval, nearly all RCs 
are well sampled (91% of the RCs in the rd sample and 
95% of those in the //opt one extend beyond 2.5 r/rd and 
0.75 r///opt, respectively) and exhibit relatively mild gra¬ 
dients. Using a range of radii rather than a single radial 
position also limits shot noise effects. The normalization 
of the RCs less extended than 3 rd (0.9 //opt) was per¬ 
formed using the available points within Ar (29 RCs in 


the rd sample were normalized using the velocity at 2.0 
rd). For each RC, we thus measured the average rota¬ 
tional velocity VAr,i] objects without a meaningful de¬ 
termination of that parameter have been excluded from 
our final sample. For a given luminosity class M, the 
template RC was obtained by scaling each RC by the 
ratio Rv^i = VAr,i/{vAr,M), where {vAr,M) is the average 
VAr,i for the class (objects with Ry^i > 2, corresponding 
to a 5cr deviation from the mean of the distribution for 
the whole sample, were discarded). The Icr dispersion 
of the VAr,i velocities is, on average, 17% of (vat.m) (it 
varies between 14% and 22% among the different lumi¬ 
nosity classes; {vAr,M) values are listed in Tables Hand 
12 presented in [Oil. We refer to the RCs scaled in this 
way as normalized RCs. 

(4) Averaging - Mean velocities and variances were com¬ 
puted (after clipping at the 2tT-level) for each grid value 
of r/rd (r///opt) for which there were at least 4 valid 
samples, and for each luminosity class. The choice of 
4 samples represents a good compromise between max¬ 
imizing the number of radial bins for the computation 
of the template RCs and keeping enough points within 
each radial bin to calculate a meaningful average veloc¬ 
ity. Since the scatter of the RCs around the template is 
artificially reduced over and near the interval Ar where 
the velocity normalization was performed, we replaced 
the velocity variance at each grid point (over the entire 
spatial range covered by the template RC, not just Ar) 
with that computed using the unnormalized RCs. In 
other words, the template RC points were computed us¬ 
ing normalized RCs, but the corresponding dispersions 
are those of the unnormalized RCs. 

(5) Smoothing - The resulting average RCs were 
Hanning-smoothed, and every other point was recorded, 
yielding a grid spacing of 0.2 rd (0.06 //opt)- As men¬ 
tioned above, each template RC point was computed by 
averaging at least 4 samples (mean velocities resulting 
from the average of 3 or less points were not recorded). 
The final template RC points were obtained by expand¬ 
ing underpopulated spatial bins (i.e., bins resulting from 
the average of 4 < A < 10 individual RC points) in or¬ 
der to include at least 10 samples. Positions and error 
bars of the expanded bins are given by the average of the 
corresponding quantities for the replaced bins. 

The binning into luminosity classes was initially done 
with a constant step AM=0.4. The bins at the faint 
end of the interval were expanded, however, to avoid 
the overlap of the corresponding template RC solutions; 
in particular, two classes Mi and Mj were combined if 
|('CAr,Mi) — ('CAr-,Mj)| < 10 km s“^. The final 10 luminos¬ 
ity classes span the /-band absolute magnitude interval 
[-24.0,-18.0] (18 objects with luminosities outside this 
range were discarded). 

The derivation of the template RCs was identical for 
the rd and Ropt parameterizations, except for the values 
of the grid spacing and the normalization interval Ar (as 
in [2] and [3] above). For the Ropt template RCs, grid 
spacing and Ar were chosen to match the corresponding 
rd ones, taking into account that the average value of 
the Ropt It’d ratio is 3.31±0.01 for the sample used to 
construct rd template RCs. 

As discussed in § ^ the inner slope of observed 
Ha/[N II] RCs depends on inclination, due to internal ex- 
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tinction in the inner regions of disks. The impact of dust 
attenuation is more severe for edge-on disks and for in¬ 
trinsically brighter systems. For this reason, 472 galaxies 
with inclination i > 80° were excluded from the analy¬ 
sis, as already mentioned. Seeing may also systematically 
bias the inner slopes of observed RCs, especially for more 
distant galaxies. The set of galaxies used to construct 
the template RCs consists of relatively nearby, large spi¬ 
rals, of median recessional velocity ~7000 km s“^ and 
isophotal radius r/ 23.5 ^ 30". We investigated a possible 
dependence of the shape of RCs on galaxy distance, and 
found no significant effect for objects in our sample. As 
another way of testing the impact of seeing on RC shape, 
we divided our sample into 3 intervals of (linear) disk 
scale length {R^ < 2 kpc, 2 kpc < i?d < 3 kpc, and 
i?d > 3 kpc) and recomputed the template RCs forcing 
the same velocity normalization obtained for the whole 
Cd sample (as done in § 0 ]for the inclination classes). We 
found that galaxies with larger i?d’s tend to have slightly 
larger rotational velocities than ones with smaller i?d’s; 
however the effect is small (well within the Icr error bars 
on the template RC velocities. In particular, the maxi¬ 
mum velocity difference is of the order of 20 km s“^ or 
less for the template RCs in the central luminosity bins, 
— 22.8 < Ml < — 20 . 8 ; there are no galaxies with i?d < 2 
kpc in the two brightest bins, nor ones with i?d > 3 kpc 
in the two faintest bins. Larger velocity differences are 
observed for the Mi = — 22.8 case, but the i?d < 2 kpc 
result is uncertain due to small number statistics) and 
limited to the r <r^ spatial interval. 

3. TEMPLATE RCS: RESULTS AND COMPARISON OF THE 
TWO PARAMETERIZATIONS 

Figure H shows the template RCs expressed as func¬ 
tions of exponential disk scale lengths rd, obtained by 
averaging the data in each of 10 luminosity intervals. 
Each curve is labeled by its mean /-band absolute mag¬ 
nitude (as listed in column 3 of Tabled described below); 
the interval Ar over which the velocities are normalized 
within each luminosity class is indicated by dotted lines. 
For each template, the outermost spatial bins have been 
expanded to include at least 10 elements; hence their 
spacing from adjacent points is not exactly 0.2 r/r^- The 
error bars are Poisson errors on the mean, calculated us¬ 
ing unnormalized RCs (see [4] in § [21). Notice that the 
error bars on the outermost points may not be accurate 
due to small number statistics. Polyex (see eq. ^ fits 
to the template RC data points are shown as solid lines. 
Our results show a smooth transition of RC shapes from 
small to more massive (luminous) systems, with the lat¬ 
ter being characterized by a steeper initial velocity rise 
and a flatter outer slope. Even in the case of the most 
luminous bins we see no convincing evidence of velocity 
decline within the spatial scales sampled here, contrar- 
ily to what expected based on the URC predictions. In 
fact, while the very outermost points of the template RCs 
show a small velocity decrease in a few cases (e.g., those 
labeled —19.37 and —22.60 in the figure), there is no sys¬ 
tematic trend that indicates a falling RC within at least 
5 disk scale lengths (6 for the central luminosity bins). 
In other words, if we were to determine the outer slope 
of each template RC by fitting a first order polynomial 
beyond 2 or 3 rd, none of the fits would have a negative 
slope. Notice also that the curve labeled —21.80 appears 


to decline beyond ~5 rd, but the two adjacent templates 
(—21.41 and —22.19) have rising outer slopes. This issue 
will be further discussed in § jS] where we compare our 
average RCs with those derived by PSS96, upon which 
the URC parameterization is built. Lastly, we point out 
that the bump at large radii seen for the lowest luminos¬ 
ity template RC is caused by a lack of faint galaxies with 
extended RCs. In fact, bright galaxies in that luminosity 
bin (i.e., those with Mi < —19.5) span all values of RC 
extent (up to 6 rd), whereas fainter ones {Mi > —19) 
have RCs that extend up to 4 rd at most (with only one 
exception, a RC with Mi = —18.6 and /?max= 4.8 rd). 

In order to illustrate the actual spread in RC shapes 
and amplitudes within each luminosity class, we include 
a version of Figure^in which the error bars represent the 
dispersion around the mean (Figure|21). Typical values of 
the standard deviation around the mean are 13 km s“^ 
for the template RCs in the lowest 4 luminosity bins, 30 
km s“^ for the one labeled —23.76, and 20 km s“^ for 
the others. 

Table n summarizes the statistics for each luminosity 
class and provides the values of the Polyex parameters 
of the template RC fits shown in Figures Q] and (21 The 
quantities listed are: 

Column (1): center of the luminosity bin. Mi, in magni¬ 
tude units. 

Column (2): width of the luminosity bin, AMi. 

Column (3): mean value of the absolute /-band magni¬ 
tude within the corresponding bin, (Mi). 

Column (4): mean velocity of the RCs, (uaf.m), mea¬ 
sured between 2rd and 3rd. This value is used to nor¬ 
malize the velocities of the individual RCs belonging to 
the same luminosity bin, as explained in § (21 (item [3]). 
Column (5): number of RCs included in each luminosity 
bin, fVAr.M- 

Columns (6)-(8): coefficients of the Polyex model (see 
eq. P) fits to the template RCs, with errors; Vq is mea¬ 
sured in km s“^, rpE is expressed in units of rd, and a 
is dimensionless. 

The bin centered at Mi= —23.8 contains only 43 objects, 
but is particularly interesting because the most luminous 
galaxies are expected to show the most declining RCs 
within a few disk scale lengths (based on URC predic¬ 
tions), implying a rather small contribution of the dark 
matter component to the internal kinematics of such sys¬ 
tems. The fastest rotators represent also important tools 
to test our un derstanding of galaxy formation and disk 
stability le.g.. lSnekkenFjl2flfl5t see also §Eli. 

The variations of the Polyex parameters with /-band 
absolute magnitude are displayed in Figure^ The trend 
seen in the top panel reflects the well-known result that 
brighter disk galaxies rotate faster (the TF relation); the 
other two panels show that RCs of brighter systems have 
steeper velocity increase within the inner regions and 
flatter outer slopes. 

The results for the template RCs expressed as func¬ 
tions of optical radii rather than disk scale lengths are 
presented in Figures 01 O 0 and Table (21 (analogous to 
Figures no El and Table n respectively). As mentioned 
in §121 the interval Ar=0.6-0.9 r//?opt for velocity nor¬ 
malization within each luminosity class corresponds to 
Ar=2-3 r/rd for an average /?opt/rd ratio of 3.3, as de¬ 
termined for the 2155 galaxies of the rd sample. The dif¬ 
ference between the corresponding values of the velocity 
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normalization (listed in column 4 of Tables ^ and 13), is 
smaller than 3 km s“^ for all the luminosity bins except 
the Ml— —23.8 one, where the discrepancy is slightly 
larger (5.4 km s“^). A quick inspection of the two ta¬ 
bles also reveals that the numbers of RCs belonging to 
the same Inminosity bin are similar, but not the same, 
for the rd and i?opt samples. Such small differences are 
caused by the exclusion of RCs that are spatially less ex¬ 
tended than 2 rd or 0.6 i?opt and by the fact that the 
.Ropt/?'d distribution is not a delta function centered on 
3.3. In fact, if we dropped the restriction on RC extent, 
or if all the RCs had i?opt=3.3 rd, the composition of 
the two samples used for template RC derivation and 
their distribution into luminosity classes would be iden¬ 
tical. In practice, for instance, a RC with a maximum 
extent i?max= 2.1 rd is part of the rd sample, but will be 
included in the i?opt sample only if i?opt/?’d< 3.5; vicev- 
ersa, a RC with i?max= 1-9 rd will be missing from the 
rd set and included in the Ropt one if i?opt/?'d< 3.2. A 
comparison between the distributions of i?opt/''’d for the 
two sets of galaxies (not shown) confirms that the i?opt 
sample has a slightly larger (smaller) fraction of objects 
with i?opt/?'d< 3.3 (> 3.3) with respect to the rd sample. 

When scaled with i?opt rather than rd, the template 
RCs have flatter outer slopes. As for the rd solutions, 
there are a few examples of templates whose outermost 
points show a small velocity decrease (Figure 0)l; how¬ 
ever, there is no clear indication of declining RCs within 
the spatial interval sampled by our data, in the sense that 
linear hts to the outer regions of the template RCs (be¬ 
yond 0.6 or 0.9 i?opt) would have flat or positive slopes. 
The only exception is the curve labeled —22.60, but again 
we notice that the adjacent templates have rising or flat 
slopes. The solutions obtained with the two parameter- 
izations are very similar, except for the tendency of the 
rd template RCs to show steeper outer slopes, as indi¬ 
cated by the larger values of the Polyex parameter a in 
the bottom panel of Figure |21 and in column (8) of Ta¬ 
ble HI compared with the corresponding results for the 
i?opt case. This effect is clearly illustrated in Figure [7| 
where the two sets of Polyex fits displayed in Figures 
n and H] are reproduced as dashed and solid lines, re¬ 
spectively. Excluding the Mi= —23.8 bin, which is also 
characterized by the poorest statistics, the template RCs 
obtained with the two parameterizations are nearly in¬ 
distinguishable within r ~1 i?opt, whereas the rd trac¬ 
ings are systematically steeper than the i?opt ones be¬ 
yond that radius. If all the RCs had the same i?opt/?’d 
ratio, the results would be identical, and in fact the rd 
templates could be obtained directly from the i?opt ones 
(or viceversa) by rescaling the horizontal axis. On the 
other hand, a distribution of Ropt/?'d values should in¬ 
troduce only scatter, and not a systematic change of RC 
slopes. In order to explain the observed bias, i?opt/?’d 
must be dependent on RC extent, in the sense that RCs 
sampled out to larger fractions of the disks must also be 
characterized, on average, by larger i?opt/?’d- Figure |H1 
shows the correlation between i?opt/^d and RC extent, 
Rmax, measured in units of disk scale lengths, for the 
2155 galaxies of the rd data set; the dotted line is at 
Ropt/'rd= 3.3, the average value for that sample, and 
the lack of points with i?max/?’d< 2 is due to a selection 
effect, as mentioned above. Our data show that for a 
fixed disk scale length, galaxies with more extended RCs 


have, on average, larger optical radii, thus explaining the 
systematic difference of slopes between the rd and i?opt 
template RCs seen in Figure [3 We investigated the ori¬ 
gin of the correlation between i?opt/?'d and RC extent 
by binning the sample into four intervals with different 
disk central surface brightness (measured from our /- 
band profiles and corrected for Galactic and internal ex¬ 
tinction, cosmological k-term, and converted to face-on 
perspective), as indicated in Figure |3 and by coding the 
objects based on their morp hological types (adopting th e 
scheme in the RC3 catalog: Fde Vaucouleurs et alJll99lD . 
Since our sample is mostly composed of Sb or S... (42%), 
Sbc (14%), and Sc (36%) galaxies, we grouped the data 
points into two classes only, Sb and earlier types (aster¬ 
isks) and Sbc and later types (open squares). Our results 
indicate that: (a) not surprisingly, the Ha emission is 
typically traced to a larger extent in the disks of galaxies 
with higher surface brightness; (b) the average value of 
Ropt / f’d increases with decreasing /tq , as the distribution 
of data points with respect to the i?opt/?'d= 3.3 line in 
each panel illustrates; (c) in the upper panel, where the 
data span the largest i?max/?'d range, the Ropt/ra ratio 
shows a clear dependence on RC extent; such trend is 
also present in the central panels, although less convinc¬ 
ingly; (d) there is no obvious dependence on morpholog¬ 
ical type, presence of an identifiable bar (not shown), or 
/-band luminosity (not shown). In conclusion, the sys¬ 
tematic difference of slopes beyond r ^i?opt between the 
two sets of template RCs in Figure [7| is largely caused by 
the high surface brightness part of our sample. In the re¬ 
gion between 4 and 6.5 i/max/j'd, where the discrepancy 
is larger, the contribution to the template RC points is 
dominated by the objects in the two upper panels of Fig¬ 
ure 13 which are characterized by an /?opt/?’d ratio that 
is systematically larger than the average of the sample, 
especially for increasing values of //max/^’d- 

In this section, we have presented template RCs ex¬ 
pressed separately as functions of exponential disk scale 
lengths or optical radii. While the description in terms 
of r<i is preferred for its more immediate interpretation 
in the context of galaxy structure and disk modeling, the 
Ropt parameterization is also provided in order to allow a 
direct comparison between PSS96’s and our results (see 
§13. From an observational point of view, individual val¬ 
ues of Cd carry signihcant errors because of the difficulty 
of fitting exponential slopes in the presence of both spi¬ 
ral structure and radially-dependent extinction. Being a 
cumulative parameter, Ropt is more straightforward to 
measure and is less subject to extinction effects. The 
complicating presence of a bulge would bias both mea¬ 
surements, in the sense that the estimated values would 
be systematically smaller than those determined for a 
galaxy with the same disk but no bulge. The impact on 
.Ropt would be direct and much larger, possibly introduc¬ 
ing morphological and/or luminosity effects. 

4. IMPACT OF INTERNAL EXTINCTION ON TEMPLATE 

RCS 

As mentioned in § 13 galaxies with inclination i > 80° 
were excluded from the sample used to derive the tem¬ 
plate RCs, in order to minimize the effects of the internal 
extinction on our results. Here we address this issue in 
more detail, and motivate our choice of the inclination 
threshold for inclusion in the final sample. 
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Several works have previously presented evidence for 
internal extinction in the inner regions of disk galaxies, 
based on the analysis of optical and near-infrared pho¬ 
tometry or spectroscopy. Extinction affects an observed 
RC by causin g its slope to appear sha llower; this effect, 
first noted bv iGoad & Robert^ ill 98 111, is illustrated by 
the radiative transfer simulations of iBosma et ;il] 1199211 
an d iBaes et al.l 1200.111 . These studies show that dust at¬ 
tenuation has a severe impact on the observed kinemat¬ 
ics of edge-on disks, even for modest optical depths, and 
that such effects are strongly reduced for galaxies that are 
more than a few degrees from edge-on, becoming entirely 
negligible for intermediate inclina tions. These results 
were c onfirmed by the analysis of iGiovanelli fc HavnesI 
1I2OOI hereafter GH02), which also showed that the 
prominence of the effect increases with disk luminosity. 
Exploiting a substantial fraction of the SFI-I--I- sample 
used in this work, GH02 showed that the kinematic scale 
length rpE of the inner slope of a RC (as obtained from 
the Polyex model in eq. ^ starts increasing noticeably 
for inclinations i > 70°, the increase being larger for 
more luminous disks. By contrast, the outer slope a of 
the observed RCs shows no dependence on inclination. 
Adopting a simple model for the dust and ga s distribu¬ 
tion in spiral disks. l\^lotto fc Giovanellil 1200411 were able 
to reproduce the observed trends in GH02 fairly well. 

In order to evaluate the impact of internal extinction 
on our results, we divided our sample, which includes 
472 highly inclined [i >80°) systems, into 3 inclination 
intervals and recomputed the template RCs. The results 
are displayed in Figure 1101 for the parameterization in 
terms of exponential disk scale lengths; only the inner re¬ 
gions are shown. Each panel corresponds to one of the I- 
band luminosity classes listed in Tabled from Mi=—23.8 
(a) to Mi=—19.0 (j); different inclination bins are indi¬ 
cated as solid (30°< i <70°), dotted (70°< i <80°), and 
dashed (80° < i <90°) lines. The velocity normalization 
of the individual RCs in each luminosity class was forced 
to be the same as that used for the whole sample, as 
listed in column 4 of Tabled Figure irni clearly shows 
that: (1) the inner slope of the most edge-on systems 
is systematically shallower, and (2) the effect is more 
pronounced for more luminous galaxies and almost dis¬ 
appears for intrinsically fainter systems. In the case of 
highly inclined galaxies, additional geometric effects can 
cause the inner slopes of observed RCs to appear shal¬ 
lower. Such effects depend on the thickness of the disk 
(higher values of disk thickness result in shallower RCs), 
which is a poorly known quantity, and can account for 
part of the difference between the highest inclination bin 
and the lower inclination ones in Figure II01 The inter¬ 
play between optical dep th at the disk center and d isk 
thickness is illustrated in iValotto fc Giovanellil 1200411 . 

Our results confirm the previous findings of GH02 (as 
expected, since they provide a different representation 
of the same effect for essentially the same data set) and 
of the other studies mentioned above, and show the 
importance of excluding the highly inclined disks from 
the derivation of the template RCs. In particular, our 
chosen threshold i = 80° allows us to effectively elimi¬ 
nate the RCs that are most affected by dust attenuation 
(and geometric effects related to the thickness of the 
disk), without discarding an excessive fraction of the 
sample. 


5. COMPARISON WITH PREVIOUS RESULTS 

A similar derivation of average RCs of disk galax¬ 
ies binned in /-band luminosity intervals was under¬ 
taken by PSS96. Their sample B included 616 RCs 
extending out to r >0.8 Rn^t (about 200 ofwhich 
beyond i?opt), obtained bv iPersic fc Sahic^ da by 
deproiecting, folding, and smoothing the Ha data of 
iMathewson. Ford, fc BiichhornI l|1 99211 . The RCs in this 
sample were normalized to the velocity measured at the 
radius encompassing 65% of the integrated light (approx¬ 
imately corresponding to 2.2 Cd or 0.7 Ropt for a pure 
exponential disk), separated into 11 luminosity intervals, 
and averaged in radial bins of size 0.1 i?opt- In order to 
characterize the velocity field of galactic disks beyond the 
optical radius, this data set was complemented by sam¬ 
ple A, a compilation of 131 RCs obtained from published 
optical or 21 cm H I-line spectroscopy. Outer gradients 
of RCs, defined as the fractional variation S of the rota¬ 
tional velocity between 1 and 2 i?opt, were calculated for 
27 RCs in sample A extending to about 2 /?opt (mostly 
H I data), and for 5 synthetic RCs obtained from the 
RCs in sample B, “suitably co-added in order to probe 
the kinematics of the outermost regions.” The variation 
of (5 as a function of R-band absolute magnitude was 
then used to extrapolate the average RCs out to 2 i?opt- 
The resulting RC tracings were fitted with a parametric 
form, obtained as the combination in quadrature of an 
exponential disk and a dark matter halo, and the fit co¬ 
efficients expressed as functions of luminosity. Based on 
this parameterization, shapes and amplitudes of observed 
RCs at any radius are completely determined by galaxy 
luminosity, implying t he existence of a univer sal RC, as 
originally proposed bv IPersic fc Sahiccil 1199111 . 

The inadequacy of the URC to model individual RCs 
has been pointed out by several authors, as mentioned 
in § n Here, we want to compare PSS96’s average RCs 
in luminosity bins (see their Figure 4) with our template 
RCs parameterized as functions of optical radii. For this 
purpose, the Polyex fits shown in Figure0]are reproduced 
in Figure El as solid lines, overlaid on the URC data 
points (kindly provided by Paolo Salucci). Since there 
seems to be a systematic offset between PSS96 and our /- 
band luminosity scales (the Mj bins for the URC tracings 
are, from top to bottom: —23.2, —22.3, —22.0, —21.6, 
—20.9, —20.5, —20.0, —19.4, —18.5), but the rotational 
velocity range sampled by our two samples is the same, 
we optimized the velocity match at the optical radius by 
shifting all t he URC data points by -1-10 km s~^. As men¬ 
tioned in Sl^ lMathewson. Ford, fc BiichhornI l|1992ll data 
set is included in our sample, but was reprocessed by us 
in the same fashion as the rest of SFI-I--I-. Therefore the 
discrepancy between absolute magnitudes might be due 
to a number of factors, including apparent magnitude 
measurements, extinction corrections, and distance esti¬ 
mates to the objects (notice that the value of the Hubble 
constant adopted by PSS96, 75 km s“^ Mpc“^, is very 
similar to the one used in this work). The comparison be¬ 
tween PSS96 and our average RCs should be carried out 
separately for the two spatial intervals below and above 
/?opt, due to the different nature of the corresponding 
URC points, as summarized at the beginning of this sec¬ 
tion. Below i?opt, where the URC points were obtained 
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by averaging RC data, the agreement is fairly good, es¬ 
pecially for the 5 lower luminosity bins (in the case of the 
curves labeled —22.2 and —22.6, one should interpolate 
between the two nearest template RCs in order to ob¬ 
tain one that matches the URC velocity at i?opt)- There 
is however a systematic tendency for the URC tracings 
to have shallower slopes, in particular for the higher lu¬ 
minosity bins (although their noise is larger), an effect 
that might be caused by the inclusion of edge-on sys¬ 
tems in PSS96 analysis (see §0. Beyond i?opt, where 
the URC points were derived from linear extrapolations 
mostly based on H I data, the agreement is good for 
the 3 lower Mi bins and becomes increasingly worse at 
higher luminosities. In particular, the URC extrapola¬ 
tions for the curves labeled Mi< —22, characterized by 
slopes that vary from marginally to strongly declining, 
are inconsistent with our Polyex fits. Interestingly, the 
results obtained by PSS96 by co-adding the RCs of their 
sample B in 5 luminosity bins show no indications of a 
velocity decrease within 2 i?opt- In particular, the curve 
labeled Mi= —22.6 in their Figure B2 (which would lie 
between the two uppermost URC tracings in Figure [mi 
has an outer gradient <5 ~ 0, indicating a flat RC slope 
between 1 and 2 i?opt- As their Figure 3 shows, the 6 
values for the 5 synthetic RCs are systematically larger 
than those derived for the individual objects in sample 
A; however, they are also represented with very large er¬ 
ror bars and said to be consistent with the others. The 
URC tracings shown in Figure I1 1 1 (with error bars re¬ 
moved for clarity) are compared with the i?opt template 
RC data points in Figure^l The URC solutions qualita¬ 
tively follow the trend of the template data in the lower 
luminosity bins, but do not recover the shapes of our av¬ 
erage RCs at higher luminosities. For instance, the data 
points of the —23.8 and —23.0 template RCs lie near the 
URC predictions, but the latter clearly do not represent 
fits to the data. In other words, both inner shapes and 
outer slopes of those template RCs, as derived from best 
fits to the data points, are not reproduced by the URC 
tracings. The disagreement is particularly evident for 
the —23.8 template RC, where a linear fit to its outer 
region (say, beyond 0.7 i?opt) is clearly inconsistent with 
the strongly declining URC slope. 

In summary, our template derivation differs from that 
of PSS96 in several details, such as the exclusion of the 
most edge-on systems, the resampling of the RCs to the 
same spatial resolution, and the normalization to a ve¬ 
locity that is both characteristic of the corresponding lu¬ 
minosity class and computed over a range of radii rather 
than at a single point. But most importantly, the larger 
size of our data set allows us to probe the velocity fields 
of the average RCs beyond 1 i?opt without reducing the 
number of luminosity bins and with no need for extrapo¬ 
lations. For radial galactocentric distances smaller than 
i?opt, where both PSS96 and our derivations are based on 
a homogeneous data set (although differently processed), 
the results are fairly consistent, when the larger scat¬ 
ter of the URC data and different selection criteria are 
taken into account. Beyond Ropt, the URC extrapola¬ 
tions qualitatively agree with our data only at low lumi¬ 
nosities, and are inconsistent for Mi< —22. PSS96 based 
their extrapolations on a linear fit to the data points in 
their Figure 3, where the outer gradient 6 is plotted as 
a function of the R-band absolute magnitude. While it 


is unclear why the outer gradients of the 5 average RCs 
built from their sample B are larger than those of the 
other points shown (mostly based on H I data), a fit to 
those 5 points would yield results in better agreement 
with our template RCs. 

We conclude that late type, high surface brightness 
spiral galaxies show no clear indication, on average, of 
declining RCs within at least 4.5 r^ (6 r^ for most of 
our Ml bins), even in the case of the most luminous sys¬ 
tems in our sample. Of course, this does not contradict 
claims of decl ining RCs at larger distances in the disks. 
Most notablv. [Ca,sertano fc va.n Corkorol (|1991jl obtained 
extended H I synthesis RCs for two highly-inclined, late 
spiral galaxies: NGC 2683 and NGC 3521. Their /-band 
absolute magnitudes, estimated from the R-band values 
using a color ind ex R —7=1.78 a ppropriate for their mor¬ 
phological type llde Jondl 199611 . and adopting Hq = 70 
km s“^ Mpc“^, are —20.6 and —22.4, respectively. Both 
galaxies exhibit declining RGs at large radii (the H I 
extends 1 5 and 12 optical ra. respectively) , and are in¬ 
cluded bv iGasertano fc van Gork(^ lll99lll in the cate¬ 
gory of “bright compact galaxies” based on their small 
values of Cd. The two RGs also show no decrease in rota¬ 
tional velocity within 6 rd and are thus consistent with 
our template results. The maximum rotation velocity of 
NGC 3521, of the order of 240 km s“^, is consistent with 
the corresponding template RC; the fact that NGC 2683, 
which is ~2 magnitudes dimmer, reaches the same max¬ 
imum velocity is somewhat puzzling but not wholly sur¬ 
prising, given the expected scatter. Extended RCs from 
H I aperture syn thesis observations were also obtained by 
iSoekken^ ll2005ll for eight of the fastest known rotators, 
with the purpose of investigating the role of disk stabil¬ 
ity in setting the observed upper limit to the rotational 
velocities of spiral galaxies. This is a very interesting is¬ 
sue, since it is unclear that the most massive disks can be 
stabilized by the same processes invoked for less massive 
syst ems; in fact, alte rnative scenarios have been proposed 
('see lSDekken^l200^ and re ferences therein), one of w hich 
predicts roughly flat RCs l|Se]1wood fc Moor^ll 999fl . For 
the eight objects in her sample (extracted from SFU—1-), 
Spekkens produced hybrid Ha-b H I RCs with extents 
ranging from 1.8 to 3.3 Ropt- All these RCs are remark¬ 
ably flat beyond Ropt (about half showing a very mild 
velocity decline for r > 1.5 — 2 Ropt) and in qualitative 
agreement with our template solutions. 

6. DISCUSSION AND CONCLUSIONS 

Taking advantage of a large sample of long-slit opti¬ 
cal RCs with available /-band photometry, we have rein¬ 
vestigated the dependence of RC shape on galaxy lu¬ 
minosity. We determined template relations by fitting 
a function, the Polyex model (eq. to the average 
RCs calculated in 10 luminosity classes, spanning the 
/-band absolute magnitude range [-24.0,-18.0]. After 
ascertaining the effect of internal extinction on our re¬ 
sults, we excluded the 472 most edge-on systems (i > 
80°) from the analysis. The template RCs are expressed 
as functions both of exponential disk scale lengths rd 
and of optical radii Ropt, the latter to allow a direct 
comparison with PSS96’s average RCs and URC predic¬ 
tions. When scaled with rd, the template relations show 
a smooth transition of RC shapes with increasing /-band 
luminosity, with the most luminous systems being char- 
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acterized by steeper inner velocity rises and flatter outer 
slopes. The i?opt templates show analogous variations of 
RC amplitudes and inner slopes with /-band luminosity, 
but their outer slopes are nearly constant. A direct com¬ 
parison of the two sets of templates shows that the results 
are very similar within /?opt; beyond that radius the rd 
templates have steeper slopes. This difference is mostly 
attributed to the high brightness galaxies, whose RCs are 
typically traced further out in the disks, and which are 
characterized by /?opt/?'d ratios larger than the average 
for the rest of the sample. We do not find convincing 
evidence for declining RCs within the spatial scales sam¬ 
pled hy our data, i.e. 4.5-6 rd, even for the most massive 
systems, in contradiction with the URC models (which, 
beyond R^t, are based on linear extrapolations). As ar¬ 
gued in there are a few examples of template RCs 
whose outermost points show a small velocity decrease, 
but there is no clear indication of declining outer slopes 
(as they would be measured from linear hts to the outer 
regions of the average RCs beyond 2 or 3 rd, or, equiv¬ 
alently, 0.6 or 0.9 i?opt)- One template RC (the one la¬ 
beled —22.60 in Figure 0 does appear to be declining 
beyond ~1.1 i?opt, but the adjacent curves have rising 
or flat outer slopes. For less luminous objects or below 
Ropt, our templates are qualitatively consistent with the 
URC tracings, when the larger scatter of PSS96’s aver¬ 
age RCs and different selection criteria (namely, the ex¬ 
clusion of edge-on galaxies from our analysis) are taken 
into account. Being signihcantly larger than the sample 
upon which the URC is based, our data set allows us 
to compute average RCs that extend beyond i?opt with- 
ont reducing the number of luminosity bins or resorting 
to extrapolation s. Our results are consistent wi th the 
declining RCs of lCa.serta.uo fc varT GorkomI 1)1991(1 , since 
the velocity decrease is observed in H I data that span 
twice the distance sampled by our optical data, and with 
the hybrid Ha-|- H I RCs of 8 fast rotators studied by 
iSDekken"^ ll2005ll . which show remarkably flat or mildly 
declining outer slopes beyond 1.5 i?opt- 
Based on our results we cannot make quantitative in¬ 
ferences about the dynamic role of dark matter in galax¬ 
ies, because optical RCs alone cannot constrain the prop¬ 
erties of dark matter halos on large scales. In fact, a de¬ 
tailed mapping of the radial distribution of visible and 
dark matter in galaxies requires the use of luminosity 
profiles and extended H I RCs in addition to optical RCs. 
Nonetheless, the comparison between Figure 0] and the 


Abadi, M. G., Navarro, J. F., Steinmetz, M., Eke, V. R. 2003, 
ApJ, 591, 499 

Avila-Reese, V., Firmani, C., &; Hernandez, X. 1998, ApJ, 505, 37 

Baes, M., et al. 2003, MNRAS, 343, 1081 

Bohm, A., et al. 2004, A&A, 420, 97 

Bosma, A., Byun, Y., Freeman, K. C., &: Athanassoula, E. 1992, 
ApJ, 400, L21 

Bullock, J. S., Dekel, A., Kolatt, T. S., Kravtsov, A. V., Klypin, 
A. A., Porciani, C., & Primack, J. R. 2001a, ApJ, 555, 240 

Bullock, J. S., Kolatt, T. S., Sigad, Y., Somerville, R. S., Kravtsov, 
A. V., Klypin, A. A., Primack, J. R., Sz Dekel, A. 2001b, 
MNRAS, 321, 559 

Casertano, S., & van Gorkom, J. H. 1991, AJ, 101, 1231 

Catinella, B. 2005, Ph.D. thesis, Cornell Univ. 

Catinella, B., Haynes, M. P., Sz Giovanelli, R., 2005, AJ, in press 
(CHG05) 

Courteau, S. 1997, AJ, 114, 2402 


URC models suggests that dark matter plays a more 
significant role in determining the internal kinematics 
within the optical disk of the most luminous systems than 
what claimed by PSS96. At high luminosity, the prohles 
of PSS96’s RCs show a strong velocity decrease between 
1 and 2 i?opt; their decomposition into a disk and a dark 
halo component in terms of the URC parameterization 
implies a modest dark matter content, whereas low lu¬ 
minosity systems are dark matter dominated. Since the 
templates RCs of the high luminosity systems do not ap¬ 
pear to decline over those spatial scales, their dark mat¬ 
ter content within the optical disk based on the URC 
predictions might be underestimated. 

The template RCs, expressed by the Polyex parameters 
as a function of galaxy luminosity, can be used to con¬ 
strain models of the circular velocity field of disks, with 
applications that range from studies of dark matter con¬ 
tent and kinematic properties of galaxies, to numerical 
simulations of disk formation and evolution within the 
current cosmological framework, and in general studies 
that, for different reasons, must rely on RC models. Our 
parameterization in terms of disk scale lengths should 
prove especially useful in the theoretical and numerical 
modeling of disk structure and evolution. 

The results presented here will be used to deter¬ 
mine statistical corrections for aperture bias of emis¬ 
sion line widths obtained with fiber spectroscopy, such 
as those that the on-going Sloan Digital Sky Survey 
(Strauss et al.i200^ is collecting for one million galaxies. 
Full rotational velocities of quiescent disk galaxies with 
apparent sizes larger than the fiber diameter can be sta¬ 
tistically recovered by simulating the impact of the finite 
aperture on the observations, where the galaxy RCs are 
modeled using our template rel ations (B. Catinella et al., 
in preparation: lCatinellall200,^ . 
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TABLE 1 

RC Distribution in Luminosity Classes and Polyex Model Fits to Template RCs Parameterized as Functions of 

Exponential Disk Scale Lengths 


Ml 

(1) 

AMi 

(2) 

( Ml ) 

(3) 

(I'Ar.M) 

(4) 

(5) 

Vb 

(6) 

The/ Td 

(7) 

a 

(8) 

-23.80 

0.40 

-23.76 

285.1 

43 

270±5 

0.37±0.02 

0.007±0.003 

-23.40 

0.40 

-23.37 

258.4 

124 

248±2 

0.40±0.01 

0.006±0.001 

-23.00 

0.40 

-22.98 

225.2 

225 

221±1 

0.48±0.01 

0.005±0.001 

-22.60 

0.40 

-22.60 

198.4 

324 

188±1 

0.48±0.01 

0.012±0.001 

-22.20 

0.40 

-22.19 

175.7 

341 

161±1 

0.52±0.01 

0.021±0.001 

-21.80 

0.40 

-21.80 

155.1 

327 

143±1 

0.64±0.01 

0.028±0.002 

-21.40 

0.40 

-21.41 

137.4 

263 

131±1 

0.73±0.02 

0.028±0.003 

-21.00 

0.40 

-21.02 

120.7 

213 

116±2 

0.81±0.02 

0.033±0.005 


0.80 

-20.48 

103.7 

203 

97±2 

0.80±0.02 

0.042±0.005 

-19.00 

2.00 

-19.37 

79.5 

92 

64±3 

0.72±0.05 

0.087±0.016 


TABLE 2 

RC Distribution in Luminosity Classes and Polyex Model Fits to Template RCs Parameterized as Functions of Optical 

Radii 


Ml 

(1) 

AMi 

(2) 

( Ml ) 

(3) 

(I’Ar.M) 

(4) 

A^Ar.M 

(5) 

Co 

(6) 

tpb/ Ropt 

(7) 

a . 

(8) 

-23.80 

0.40 

-23.76 

290.5 

44 

275±6 

0.126±0.007 

0.008±0.003 

-23.40 

0.40 

-23.37 

260.6 

130 

255±2 

0.132±0.003 

0.002±0.001 

-23.00 

0.40 

-22.98 

227.8 

226 

225±1 

0.149±0.003 

0.003±0.001 

-22.60 
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Fig. 1. — Template RCs parameterized as functions of exponential disk scale lengths. Each curve is labeled on the right by its mean 
/-band absolute magnitude (Hq = 70 km s~^ Mpc“^). The final sample includes 2155 RCs extending beyond 2 and with inclination 
to the line of sight i <80°. The vertical, dotted lines show the interval over which the velocity normalization was performed (see 
error bars are Poissonian errors on the mean. Polyex fits to the data points are indicated by solid lines; the fit coefficients are presented in 
Figure [3 and Table [Tl 
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:e RCs shown in Figure^ are displayed here in separate panels, with error bars representing the dispersion around 
3S indicate the Polyex fits from Figure [Tl 
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Fig. 3.— Polyex model coefficients of the fits shown in Figure^ plotted as functions of J-band absolute magnitude. See also Table [Tl 
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Fig. 4. — Same as Figure [Tl for template RCs parameterized as functions of optical radii. This sample includes 2169 RCs extending 
beyond 0.6 Ropt, and with inclination i <80°. Polyex model coefficients are presented in Figure!^ and Table |21 
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irel^for the -Ropt parameterization. The solid lines indicate the Polyex fits from Figure|^ 
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Fig. 6. — Polyex model coefficients of the fits shown in Figure|4| see also Tabled The vertical scales are the same as those in Figure El 
except for the central panel. 
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f Ropt (solid lines) or (dashed). The upper horizontal scale 
e of the Ropt/^d ratio for the galaxies of the sample. 
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Fig. 8.— Correlation between optical radius and RC extent, both expressed in units of J-band exponential disk scale lengths, for the 
sample. The dotted line is at Ropt/?’d= 3.3, the average value for this data set. The trend observed here accounts for the systematic 
difference of outer slope between the two sets of template RC shown in FigurelTI 





6 

4 

2 

6 

4 

2 

6 

4 

2 

6 

4 

2 

le 


Catinella et al. 




18.5 < ^19.5 



^ /Ltg > 19.5 
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Figure[3 with the sample divided into four bins of central surface brightness, /io, as indicated in each panel. Galaxies 
:ype Sb or earlier are represented by asterisks, those classified as Sbc or later types are marked as open squares. 
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FlCs divided into 3 inclination intervals: 30°< i <70°(solid lines), 70°< i <80°(dotted), and 80°< i <90°(dashed), 
orrespond to the luminosity classes Mi=—23.8 to —19.0 in Table [Tl Only the inner regions of the RCs are shown. 
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Fig. 11. — Polyex fits (solid lines) to the template RCs shown in Figure|^are here compared to the average RCs obtained by PSS96 and 
kindly provided by Paolo Salucci. All PSS96 data points have been shifted by +10 km s“^ to match the velocities of the template curves 
at Ropt • 
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Fig. 12. — The URC tracings in Figure FTTl are here reproduced as solid lines, superimposed on the -Ropt template RC data from Figurel^ 
















